Two-Step Model of Fusion for Synthesis of Superheavy Elements 



(N 
O 

o 

(N 
o 

Q 

o 



(N 
> 

(N 
O 
(N 

(N 
O 

i 

=5 



X 



Caiwan Shen 1 ' 2,4,6 , Grigori Kosenko 3 ' 4 and Yasuhisa Abe 5 
1 China Institute of Atomic Energy, P.O.Box 275(18), Beijing 102413, China 
2 INFN-LNS, U Via S. Sofia, 1-95123 Catania, Italy 

3 Department of Physics, University of Omsk, Omsk, Russia 

4 RIKEN, 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan 

5 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, Japan 
6 Center of Theoretical Nuclear Physics, National Laboratory of Lanzhou, Lanzhou, China 

A new model is proposed for fusion mechanisms of massive nuclear systems where so-called fusion 
hindrance exists. The model describes two-body collision processes in an approaching phase and 
shape evolutions of an amalgamated system into the compound nucleus formation. It is applied to 
48 Ca- induced reactions and is found to reproduce the experimental fusion cross sections extremely 
well, without any free parameter. Combined with the statistical decay theory, residue cross sections 
for the superheavy elements can be readily calculated. Examples are given. 
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How many elements exist in the nature or what is the 
heaviest element has been an intriguing question since 
the periodic table was proposed for the chemical ele- 
ments. The heaviest element that exists in the nature is 
now known to be Uranium with atomic number Z being 
92. But the discovery of the magic numbers in atomic 
nuclei and their understanding by the shells of nucle- 
onic motion[l] suggest that much heavier atomic nuclei 
might exist, stabilized by extra-bindings due to possible 
shells next to the largests known, i.e., Z=82 and N=126. 
Actually, many theoretical calculations have been made, 
predicting the next double closed shell nucleus to be 
with Z=114, 120, or 126 and N=184[2]. Naturally, enor- 
mous experimental efforts have been devoted to finding 
out traces of existence of the corresponding superheavy 
atomic nuclei and to synthesizing them with nuclear reac- 
tions, especially with heavy-ion fusion reactions [3] . But 
what combination of ions is favorable as entrance chan- 
nels and what incident energy is the optimum for residues 
are not predicted well, and thus, the experiments have 
been performed according to the results of systematic 
studies done so far. This is due to the lack of our knowl- 
edge of reaction mechanisms. 

Based on the theory of compound nucleus reactions, 
the residue cross sections are given as follows, 

a rcB = 7rX 2 Ej(2 J + 1) • PjU n (E c . m .) • P S J UIV (E*), (1) 

where X is the inverse of the wave number and J is the 
total angular momentum quantum number. Pf US ion and 
^surv denote the fusion and the survival probabilities, 
respectively. The latter is given by the statistical the- 
ory of decay, i.e., by competitions between neutron emis- 
sion and fission decay. Essentially unknown is the fusion 
probability, i.e., fusion mechanism of massive systems, 
although there are ambiguities in the parameters in the 
properties of heavy and superheavy nuclei which give rise 
to uncertainties in calculating the survival probability. 

In lighter systems, the fusion probability is well deter- 
mined by the barrier defined with the Coulomb and the 



nuclear attraction between nuclei in the entrance chan- 
nel, but in massive systems, the situation is not so sim- 
ple. It has been well known experimentally that there is 
the fusion- hindrance [4] , which is often described with so- 
called extra-push energy which is required for a system to 
fuse in addition to the barrier height [5]. A physical ori- 
gin or mechanism is not yet well clarified. There are two 
possible interpretations proposed. They both attribute it 
to energy dissipations; one is due to the dissipation of the 
initial kinetic energy during two-body collisions passing 
over the barrier [6], while the other is due to the dissi- 
pation of the energy of collective motions which would 
lead an amalgamated system to the spherical compound 
nucleus [5]. It is natural to consider that both mech- 
anisms exit, though we don't know a priori which one 
dominates in which situation. We, thus, propose a new 
theoretical framework for fusion, i.e., a two-step model 
which incorporates both of them properly [7] . 

In the approaching phase of passing over the Coulomb 
barrier, we describe the system as collision processes un- 
der frictional forces, up to the contact point of two in- 
cident nuclear matters and then, we describe dynamical 
evolutions of the amalgamated mono-nuclear system to- 
ward the spherical shape under frictional forces acting 
in collective motions of excited nuclei. As is given be- 
low, both dynamical processes are described by Langevin 
equations which include random forces associated to the 
respective frictions. It would be worth to mention here 
that the fluctuations due to the random forces are cru- 
cially important in problems of small probability such 
as in syntheses of the superheavy elements (SHE), be- 
cause we have to investigate cases where mean trajecto- 
ries never reach the spherical shape. Another point to 
be mentioned is that since the two steps are connected 
successively, the results of the first step not only gives a 
probability for incident ions to stick to each other (stick- 
ing probability P s tick) but also gives initial conditions for 
the second step. Thus, the method of the connection 
from the first to the second steps is natural, which is nei- 
ther related to the diabaticity nor the adiabaticity. It is 
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completely new and could be called as "statistical" , as 
will be seen below. In massive systems, there is a con- 
ditional saddle point, or a ridge line between the amal- 
gamated configuration and the spherical shape on the 
potential energy surface calculated with the liquid drop 
model (LDM), which could be considered to be another 
barrier inside and makes most trajectories to return back 
to re-separation (quasi-fission, etc.), i.e., gives rise to a 
small probability for forming the spherical shape (forma- 
tion probability Pf or m)- Thus, the fusion probability is 
given by the product of the two probabilities, 



^"fusion C^c.m.) 



-^stick (-^c.m.) ' -fform^c.r 



(2) 



In order to realize the model, we employ the surface fric- 
tion model (SFM)[8] for the approaching phase and the 
one-body wall-and-window formula[9] of the dissipation 
for the shape evolutions, i.e., for the second step. 

As for the approaching phase, the equation of motion 
is only for the radial degree of freedom and the orbital 
angular momentum, and is given below, 

dr 1 
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(L(t) - L at ) + R T {t), 



(3) 



where \i is the reduced mass of the collision system, and 
V is the sum of the Coulomb potential V c and the nuclear 
potential V n . Cj(r) is the radial or the tangential friction 
coefficient which is assumed to have the following form 
factor, 



Ci(r) = K° ■ (dV n /drY 



(4) 



where K° = 0.035 and K% = 0.0001 in unit of 
10~ 21 s/MeV. Ri denotes a random force associated with 
the friction for i = r (radial) or T (tangential), and as- 
sumed to be Gaussian, and to satisfy the following prop- 
erty, 

(Ri(t)) = 0, 

(Ri(t) ■ Rjit')) = 25«<J(t - t') ■ d (r(t))T J (t), (5) 



where the last equation is the fluctuation-dissipation the- 
orem with temperature T J (t), J being equal to a total 
angular momentum of the system, i.e., an incident or- 
bital angular momentum L. L st denotes the limiting or- 
bital angular momentum under the friction, which is so- 
called the sliding limit in the SFM and is equal to 5/7- L. 
We calculate many trajectories over relevant impact pa- 
rameters and obtain probabilities for their reaching the 
contact point, respectively. Fig. la shows the calculated 
sticking probability for L=0 for the case of 48 Ca- 238 U 
system. Incident energy is given relative to the barrier 
height. It is readily seen that at energies just above the 
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FIG. 1: Results on 48 Ca- 238 U system with SFM. The stick- 
ing probability for L=0 is shown in (a) (curve is for an eye- 
guide), the radial momentum distribution in (b) in unit of 
l(F 21 sec-MeV/fm and the average orbital angular momentum 
v.s. the relative distance in (c). 



barrier there is almost no probability. This is due to the 
fact that the form factor assumed in SFM stretches over 
outside the barrier top position in massive systems. The 
results already appear to explain the fusion hindrance 
and at least partially the extra-push energy, while the 
second step is also expected to give rise to an additional 
contribution. In order to know the physical situation at 
the contact point, we analyze the radial momentum dis- 
tribution as well as that of the orbital angular momen- 
tum. The radial momentum distribution is found to be 
almost purely Gaussian, as shown in Fig. lb. Its width is 
in consistence with the temperature of the heat bath of 
nucleons which is supposed to absorb the initial kinetic 
energy through the friction force. The example shown is 
for L = 0, but the other angular momentum cases behave 
in the same way. Therefore, the calculated distribution 
S J (po, £7 c .m.) can be expressed as follows for each angular 
momentum, 



S J (po,£c.m.) 



P s i ck (E c . m .)-g J (p ,p J ,T J ), 



(6) 



where the normalized Gaussian distribution 
9 J {Po, Po> Tq) is given generally so as to include 
an average mean momentum left (Pq) which is almost 
equal to zero in the present case. This distribution is 
used as the initial inputs to the dynamical evolutions 
in the second step, i.e., to Eq. (8) below. T J denotes 
the temperature of the amalgamated system. The total 
energy available for the compound nucleus E* is written 
by the energy conservation as follows, 



E* = £ c . m . + Q = Vq — -Eshoii + £o + ko, 



(7) 



where Q denotes the Q- value of the fusion reaction. E s h e n 
the shell correction energy of the ground state, Vo the 
LDM potential energy of the contact point, So the in- 
trinsic excitation, and ko the radial kinetic energy left 
at the contact point. The latter two are averagely given 
as ao • Tq 7 and (pg ) 2 /2yU + \Tq , respectively with the 
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level density parameter ao which is calculated according 
to Toke and Swiatecki [10]. The orbital angular momen- 
tum is also analyzed. The average value is plotted as 
a function of the radial distance in Fig.lc. It is seen 
that it approaches to the dissipation limit L st about the 
contact point, which indicates that the incident system 
reaches the sticking limit, if the rolling friction is prop- 
erly taken into account. We, thus, can consider that the 
relative motion is completely damped and reaches the 
thermal equilibrium with the heat bath at the contact 
point, i.e., that the incident ions form an amalgamated 
mono-nuclear system, the probability of which depends 
on the incident energy and is extremely small just above 
the barrier. It should be noticed here that p'l = does 
not always hold, for example not in 100 Mo- 100 Mo system 
etc. 

Subsequent shape evolutions of the pear-shaped mono- 
nucleus formed with the incident ions are described by 
the multi-dimensional Langevin equation which is the 
same as that used for dynamical studies of fission[ll], 
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FIG. 2: Examples of the trajectories are displayed with the 
same initial radial momentum being equal to zero. Random 
force gives rise to a variety of the trajectories. The circle in the 
upper right corner corresponds to the touching configuration 
reached by the first step, from which dynamical evolutions 
of shape start. (Ro being the radius of the spherical ground 
state) 

where summation is implicitly assumed over repeated suf- 
fixes. The collective mass tensor my is the hydrodynam- 
ical one and the potential U J is calculated by the finite 
range LDM with two-center parameterization of nuclear 
shapes [12], added with the rotational energy of the sys- 
tem calculated with the rigid body moment of inertia. 
The random force Ri{t) is again gaussian with the nor- 
malization 2, and the tensor g^j is now related to the 
friction tensor jij, as is given in the last equation, i.e., 



the generalized Fluctuation-Dissipation theorem in the 
multidimensional case. The friction tensor is calculated 
with the wall-and- window formula [9]. The temperature 
T J of the heat bath is better to be taken to be that at 
the conditional saddle point, but is approximated with 
that at the contact point, i.e. Tq . They are close to 
each other for the 48 Ca- induced reactions. In the present 
calculations we only use the relative distance R and the 
mass asymmetry coordinate a with the other degrees of 
freedom being frozen. For example, the neck parameter 
is taken to be 0.8, based on our experiences that it does 
not change so much during passing over the conditional 
saddle point in the three-dimensional calculations. Fig. 2 
shows examples of the trajectories on the LDM potential 
for 48 Ca- 238 U system for initial radial momenta and thus 
initial energies being equal to zero. Calculations of many 
trajectories, starting with various initial radial momenta 
give a distribution of formation probability F J (p ,T J ). 
By making a convolution of it with the Gaussion distri- 
bution of the initial momentum g J (po,Pq , Tq), we obtain 
the formation probability Pf or m 

Pform(£c.m.) = J dp F J (p ,T J ) ■ g J (p ,P^ Tq). (9) 

Fig. 3a shows the calculated formation probability for 
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FIG. 3: Calculated formation and fusion probabilities are 
shown in (a) and (b), respectively. 

48 Ca- 238 U system for L=0 and 30. In Fig.3b, the final 
fusion probability is plotted versus incident energy. At 
the first glance, the decreasing energy dependences seem 
to be peculiar, but the energy dependence of the passing- 
over probability under friction delicately depends on the 
strength of friction and the incident momentum. Ac- 
tually, slightly weaker friction gives rise to an increas- 
ing energy dependence. A detailed analysis with the 1- 
dimensional model will be given elsewhere [15]. It should 
be also mentioned here that the present model is com- 
pletely classical, and thus there is no quantum tunnel- 
ing effect included, which limits the lowest energy to be 
reached. 

Fusion cross sections are calculated with the fu- 
sion probability as usual, (Tf us i on = ttX Sj(2J + 1) • 



4 




180 200 220 



FIG. 4: Calculated excitation functions of fusion reactions 
for 4S Ca- 244 Pu, - 24S Cm and - 252 Cf systems, together with the 
available experimental data [13] for GSI and [14] for Dubna. 

TABLE I: Calculated maximum residue cross sections of the 
three systems 48 Ca- 244 Pu, - 248 Cm and - 252 Cf are summarized. 
The factor 1/3 for Moller masses is rather arbitrarily chosen. 
(o- max :pb, E*\ MeV) 
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244 Pu 


Liran —0.23 


0.018 


30.6 


0.018 


36.5 
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7.39 


30.1 


6.00 


35.3 




Experiment 


« 1 


Blab =236 


w 1 


£lab = 236 


248 Cm 


Liran —1.37 


0.254 


31.1 


0.045 


37.8 
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4.56 
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2.98 


35.6 
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0.6 


35.8 


2b2 Cf 


Liran -3.24 


1.057 


32.7 


0.095 


38.2 




M611er/3 -2.41 


0.216 


28.8 


0.086 


33.5 



-ffusionC^'c.m.), and are shown in Fig. 4 for the four systems 
with Ca beam, together with some measured cross sec- 
tions [14]. It is extremely surprising that the calculations 
well reproduce the experiments without any adjustment 
of the model parameters. Experimental measurements 
are highly desirable in other heavy systems for compar- 
isons with the present calculations. 

In order to show that we are ready for calculations 
of residue cross sections for SHE, we give examples for 
Z=114, 116 and 118, by the use of P SU rv calculated with 
HIVAP [16]. Actually, the shell correction energies are 
the most crucial quantities in residue calculations, be- 
cause they effectively give the fission barriers for SHEs. 
And they are not yet firmly predicted, and thus we take 
those by Moller and Liran [17] as typical examples of 
mass predictions, and compare with the recent Dubna 
experiments [18], which are given in Table I. 

In brief, the new two-step model has been found to 
be extremely successful in reproducing the available fu- 
sion data of 48 Ca induced reactions. By combining the 
present fusion probabilities with the standard statistical 
decay calculations, we have obtained residue cross sec- 
tions for Z=114, 116, and 118, which are in a reasonable 
agreement with the recent Dubna experiments, but with 
rather small shell correction energies, much smaller than 
previously thought. A systematic study of residue cross 
sections are being made. Furthermore, the model is now 
being applied to other massive systems, such as 100 Mo- 
100 Mo etc. 
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